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Abstract 

Adaptive approximation (or interpolation) takes into account local variations in the behavior of 
the given function, adjusts the approximant depending on it, and hence yields the smaller error of 
approximation. The question of constructing optimal approximating spline for each function proved 
to be very hard. In fact, no polynomial time algorithm of adaptive spline approximation can be 
designed and no exact formula for the optimal error of approximation can be given. Therefore, 
the next natural question would be to study the asymptotic behavior of the error and construct 
asymptotically optimal sequences of partitions. 

In this paper we provide sharp asymptotic estimates for the error of interpolation by splines on 
block partitions in IR'*. Wc consider various projection operators to define the interpolant and provide 
the analysis of the exact constant in the asymptotics as well as its explicit form in certain cases. 
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1 Introduction 



The goal of this paper is to study the adaptive approximation by interpolating splines defined over block 
partitions in R'^. With the help of introduced projection operator we shall handle the general case, 
and then apply the obtained estimates to several different interpolating schemes most commonly used in 
practice. 

Our approach is to introduce the "error function" which reflects the interaction of approximation 
procedure with polynomials. Throughout the paper we shall study the asymptotic behavior of the ap- 
proximation error and, whenever possible, the explicit form of the error function which plays a major 
role in finding the constants in the formulae for exact asymptotics. 



1.1 The projection operator 

Let us first introduce the definitions that will be necessary to state the main problem and the results of 
this paper. 

We consider a fixed integer d> 1 and we denote by a; = (xi, • • • , Xa) the elements of M''. A block R 
is a subset of M** of the form 

l<i<d 

where Ui < bi, for all 1 < i < d. For any block i? C M'', by Lp{R), 1 < p < oo, we denote the space of 
measurable functions / : i? — )• R for which the value 



\f{x)\Pdxj , if l<p<oo, 

esssup{|/(a;)| : x G R}, if p = oo. 



is finite. We also consider the space C^{R) of continuous functions on R equipped with the uniform norm 
II • ||Lo5,(fi)- We shall make a frequent use of the canonical block I'', where I is the interval 



I := 



1 1 
2' 2 



Next we define the space V := C°{I'^) and the norm || • jj^ := || • ||z,^(id). Throughout this paper we 
consider a linear and bounded (hence, continuous) operator I : V ^ V. This implies that there exists a 
constant Ci such that 

\\ I u\\v <Ci\\u\\v ioT all u€V. (1.1) 
We assume furthermore that I is a projector, which means that it satisfies 

IoI = I. (1.2) 

Let R be an arbitrary block. It is easy to show that there exists a unique a;o € M."^ and a unique diagonal 
matrix D with positive diagonal coefficients such that the transformation 

(j){x) := Xo + Dx satisfies 0(1'') = R. (1.3) 

The volume of R, denoted by |i?|, is equal to det(D). For any function / e C^{R) we then define 

(1-4) 

Note that 

11/ - Ir /IU,(«) = (det D)T.\\fo^-l{fo 4>)\\L,(l^y (1.5) 

A block partition 7?. of a block Rq is a finite collection of blocks such that their union covers Rq and 
which pairwise intersections have zero Lebesgue measure. If 7?. is a block partition of a block Rq and 
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if / G C'~'{Ro), by Itj / e ioo(-Ro) we denote the (possibly discontinuous) function which coincides with 
Ifl / on the interior of each block R E TZ. 

Main Question. The purpose of this paper is to understand the asymptotic behavior of the quantity 

11/ - Ikjv /IkpCflo) 

for each given function f on Rq from some class of smoothness, where (TZn)n>i is a sequence of block 
partitions of Rq that are optimally adapted to /. 

Note that the exact value of this error can be explicitly computed only in trivial cases. Therefore, the 
natural question is to study the asymptotic behavior of the error function, i.e. the behavior of the error 
as the number of elements of the partition TZ^ tends to infinity. 



Most of our results hold with only assumptions (1.1 1 of continuity of the operator I, the projection 



axiom (1.2), and the definition of 1^ given by (1.4 1. Our analysis therefore applies to various projection 
operators I, such as the L2-orthogonal projection on a space of polynomials, or spline interpolating 
schemes described in 91.41 



1.2 History 

The main problem formulated above is interesting for functions of arbitrary smoothness as well as for 
various classes of splines (for instance, for splines of higher order, interpolating splines, best approximating 
splines, etc.). In the univariate case general questions of this type have been investigated by many authors. 
The results are more or less complete and have numerous applications (see, for example, [12] )• 

Fewer results are known in the multivariate case. Most of them are for the case of approximation by 
sphnes on triangulations (for review of existing results see, for instance [HI El [2 [3 [T^] ) . However, in 
applications where preferred directions exist, box partitions are sometimes more convenient and efficient. 

The first result on the error of interpolation on rectangular partitions by bivariate splines linear in each 
variable (or bilinear) is due to D'Azevedo [8 who obtained local (on a single rectangle) error estimates. 
In [3] Babenko obtained the exact asymptotics for the error (in Li, L2, and Loo norms) of interpolation 
of C^(I'') functions by bilinear splines. 

In [5^ Babenko generalized the result to interpolation and quasiinterpolation of a function / e C'^{1'^) 
with arbitrary but fixed throughout the domain signature (number of positive and negative second-order 
partial derivatives). However, the norm used to measure the error of approximation was uniform. 

In this paper we use a different, more abstract, approach which allows us to obtain the exact asymp- 
totics of the error in a more general framework which can be applied to many particular interpolation 
schemes by an appropriate choice of the interpolation operator. In general, the constant in the asymp- 
totics is implicit. However, imposing additional assumptions on the interpolation operator allows us to 
compute the constant explicitly. 

The paper is organized as follows. Section [l . 5| contains the statements of main approximation results. 
The closer study of the error function, as well as its explicit formulas under some restrictions, can be 
found in Section [2] The proofs of the theorems about asymptotic behavior of the error are contained in 
Section [Sj 



1.3 Polynomials and the error function 

In order to obtain the asymptotic error estimates we need to study the interaction of the projection 
operator I with polynomials. 

The notation a always refers to a d-vector of non-negative integers 

a = (ai, • • • ,ad) e K^J.. 

For each a we define the following quantities 



:— Q!;, a! :— a^!, max(a) :— max 

We also define the monomial 



l<i<d l<i<a 



l<i<d 
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where the variable is X = (Xi, Xd) G M^. Finally, for each integer fc > we define the following three 
vector spaces of polynomials 



TPk 

n 



= Vect{X" 
= Vectjx" 
= Vectjx" 



|a| < k}, 

max(a) < k and jaj < + 1}, 
max(a) < fc}. 



(1.6) 



Note that clearly dini(P^* 



Furthermore, 



{k + l^. In addition, a classical combinatorial argument shows that 

, ^ , k + d 
dimPfc = ( ^ 



dimP^. = dimPfc+i — d = 



k + d+1 
d 



-d. 



By Vi we denote the image of I, which is a subspace of = C^{¥'-). Since I is a projector (1.2), we 
have 

= {!(/): f^V}^{f^V: / = !(/)}. (1.7) 
From this point on, the integer k is fixed and defined as follows 

k = fc(I) max{fc' > : Pfc- C Vi] (1.8) 

Hence, the operator I reproduces polynomials of total degree less or equal than k. (If fc = oo then we 

obtain, using the density of polynomials in V and the continuity of I, that I(/) / for all f E V. We 
exclude this case from now on.) 

In what follows, by m we denote the integer defined by 



m = rn(I) := fc + 1, 



(1.9) 



where k — fc(I) is defined in ( 1.8 1. By we denote the space of homogeneous polynomials of degree m 

Mm := Vect{X" : |a| = m}. 

We now introduce a function Kj on Hm , further referred to as the "error function" . 

Definition 1.1 (Error Function) For all tt e Mm 

Kiiir):^ M Wn-lnnU^^j^j, (1.10) 
I j'ti i 

where the infimum is taken over all blocks R of unit d-dimensional volume. 

The error function K plays a major role in our asymptotical error estimates developed in the next 
subsection. Hence, we dedicate §2 to its close study, and we provide its explicit form in various cases. 

The optimization ( 1.10 1 among blocks can be rephrased into an optimization among diagonal matrices. 
Indeed, if \R\ = 1, then there exists a unique xq G K'' and a unique diagonal matrix with positive 
coefficients such that R = (f>{l'^) with (p{x) — Xq + Dx. Furthermore, the homogeneous component of 
degree m is the same in both tt o cjj and tt o D, hence tt o cj) — n o D € P^ (recal that m = k + 1) and 
therefore this polynomial is reproduced by the projection operator I. Using the linearity of I, we obtain 

IT o <j) — I(7r o (j)) = TT o D — 1{tt o D) . 

Combining this with (1.5), we obtain that 

Klin) 



inf |l7roi:)-I(7roi:))||i (j<i), 

dctZ)— 1 ^ 
D>0 



(1.11) 



where the infimum is taken over the set of diagonal matrices with non-negative entries and unit determi- 
nant. 
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1.4 Examples of projection operators 

In this section we define several possible choices for the projection operator I which are consistent with 
(1.8 1 and, in our opinion, are most useful for practical purposes. However, many other possibilities could 
be considered. 

Definition 1.2 (^2(1'^) ortiiogonal projection) We may define !(/) as the L2{1'^) orthogonal projec- 
tion of f onto one of the spaces of polynomials P^, IP*/, or IPl* defined in (1.6). 

If the projection operator I is chosen as in Definition [L2] then a simple change of variables shows that 
for any block i?, the operator I/? defined by (1.4) is the L2{R) orthogonal projection onto the same space 
of polynomials. 

To introduce several possible interpolation schemes for which we obtain the estimates using our 
approach, we consider a set C/fc C I of cardinality #(C/fc) = A: + 1 (special cases are given below). For any 
u = (ui, • • • Ud) e we define an element of F^* as follows 



l<i<d 



n 



Clearly, y.^{u) = Aiu(wi,-- - ,Ud) = I and Aiu(v) = /iu(wi,--- ,1;^) if v = {vi,--- ,Vd) € and 
V 7^ u. 

It follows that the elements of B := (/iu)ueL"* linearly independent. Since #(i?) = = 
(fc + l)'^ = dim(Fr), 5 is a basis of TPl* . 

Therefore, any element of /i G can be written in the form 

fi{X) = ^uA^u(^). 

It follows that there is a unique element of /i e such that /i(u) = /(u) for all u G U^. We define 
I / := /i, namely 

(I/)(X):= /(uK(x)ePr. 

We may take Uk to be the set of A; + 1 equi-spaced points on I 

1 



Uk = <- 



n 
k 



<n< k 



(1.12) 



We obtain a different, but equally relevant, operator I by choosing Uk to be the set of Tchebychev points 
on I 



1 



cos 



UTT 

T 



< n < fc 



(1.13) 



Different interpolation procedures can be used to construct I. Another convenient interpolation scheme 
is to take 

and I(/) = / on a subset of U^. This subset contains dimPjl points, which are convenient to choose first 
on the boundary of I'* and then (if needed) at some interior lattice points. Note that since dimP^ < 
#{U^) = (fc + 1)'', it is always possible to construct such an operator. 

If the projection operator I is chosen as described above, then for any block R and any / e C^{R), 
^r{I) is the unique element of respective space of polynomials which coincides with / at the image 4'{p) 
of the points p mentioned in the definition of I, by the transformation described in ( 1.3 1. 
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1.5 Main results 



In order to obtain the approximation results we often impose a slight technical restriction (which can be 
removed, see for instance 2J) on sequences of block partitions, which is defined as follows. 

Definition 1.3 (admissibility) We say that a sequence {'JZn)n>i of block partitions of a block Rq is 
admissible if #{TZn) < N for all N >1, and 



sup ( sup diam(i?) ) < oo 
N>i V ReiiN J 



(1.14) 



We recall that the approximation error is measured in Lp norm, where the exponent p is fixed and 
1 < p < oo. We define r G (0, oo) by 

1 m 1 

(1.15) 



1 TO 1 

rap 

In the following estimates we identified d™f{x) with an element of Hm according to 



d^'fix) 



E 



a"7(a;) X" 
dx" al 



(1.16) 



We now state the asymptotically sharp lower bound for the approximation error of a function / on an 
admissible sequence of block partitions. 

Theorem 1.4 Let Rq be a block and let f G C™(i?o)- For any admissible sequence of block partitions 
{TZn)n>i of Ro 

'd^r 



liminfiVrf ||/-l7j„/|U^(fl„) 



> 



Kr 



The next theorem provides an upper bound for the projection error of a function / when an optimal 
sequence of block partitions is used. It confirms the sharpness of the previous theorem. 

Theorem 1.5 Let Rq be a block and let f G C""(i?o)- Then there exists a (perhaps non-admissible) 
sequence {TZn)n>i, it^Ti-N < N, of block partitions of Rq satisfying 



limsupiV" 11/ - Ik„ /|Up(i?o) < 

AT— >-oo 



K, 



d'^f 



(1.17) 



Furthermore, for all s > there exists an admissible sequence {TI%)n>i of block partitions of Rq 
satisfying 

limsup7V^|l/-lK^^/|U^(H„)< Kii^) +e. (1.18) 

An important feature of these estimates is the "lim sup" . Recall that the upper limit of a sequence 
{un)n>No is defined by 

limsupiiAT := lim sup m„, 

JV-foo N^oOn>N 

and is in general strictly smaller than the supremum sup^>jV(, un- It is still an open question to find an 
appropriate upper estimate of supjY>Ar^ A^^||/ — Ik„ (Aq) when optimally adapted block partitions 
are used. 

In order to have more control of the quality of approximation on various parts of the domain we 
introduce a positive weight function il G C'^{Rq). For 1 < p < oo and for any u G Lp(i?o) as usual we 
define 

\\u\\L^{Ro,n) ■■= \\un\\L^(^R^y (1.19) 



Remark 1.6 Theorems 



1.4 



1.5 



and 



1.7 below also hold when the norm 



, Wlj rq) 



(resp 



replaced with the weighted norm \\ ■ \\Lp(Ra,n) (resp \\ ■ \\l^(Rq,q.)) defined in {1.19). 



\L^{Ro 



)) is 
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In the following section we shall use some restrictive hypotheses on the interpolation operator in order 
to obtain an explicit formula for the shape function. In particular, Propositions |2.7[ |2.8[ and equation 
(2.20) show that, under some assumptions, there exists a constant C = C(I) > such that 



1 ^ fd-^f 



< 



\ 



n 



l<i<d 



dx]' 



< CKi 



7I 



These restrictive hypotheses also allow to improve slightly the estimate (|1.18|) as follows 



Theorem 1.7 // the hypotheses of Proposition 



2.1 or 



2. 8 hold, and if Kj 



> everywhere on 



Rq, then there exists an admissible sequence of partitions {TIn)n>i which satisfies the optimal estimate 



1.17). 



The proofs of the Theorems 1.4 1.5 and 1.7 are given in J|3j Each of these proofs can be adapted to 
weighted norms, hence establishing Remark 1.6 Some details on how to adapt proofs for the case of 
weighted norms are provided at the end of each proof. 



2 Study of the error function 

In this section we perform a close study of the error function Kj, since it plays a major role in our 
asymptotic error estimates. In the first subsection §2.1| we investigate general properties which are valid 
for any continuous projection operator I. However, we are not able to obtain an explicit form of Kj 



under such general assumptions. Recall that in Q.4 we presented several possible choices of projection 



operators I that seem more likely to be used in practice. In §2.2| we identify four important properties 
shared by these examples. These properties are used in |2.3| to obtain an explicit form of Kj. 

2.1 General properties 

The error function K obeys the following important invariance property with respect to diagonal changes 
of coordinates. 

Proposition 2.1 For all tt G IHm and all diagonal matrices D with non-negative coefficients 

Ki (tt o D) = (dot D)'^Ki{ti). 

Proof: We first assume that the diagonal matrix D has positive diagonal coefficients. Let Z? be a 
diagonal matrix with positive diagonal coefficient and which satisfies detZ? = 1. Let also tt G Hm. Then 

TTo{DD) = TTo ((detD)3i)) = {dctD)'^Tro D, 

where D := {det D)^i DD satisfies detl) = det D = 1 and is uniquely determined by D. According to 
( |1.11[ ) we therefore have 

KiiiroD) = inf \\7r o (DD) - ^rr o {DD))\\ ^^^j.) 

dct^D=l 
D>Q 

= (detT^)? inf ||7roD-I(^o^)||^ (I.) 

detLI = l 

= [dei D)^ Ki{-k), 

which concludes the proof in the case where D has positive diagonal coefficients. 

Let us now assume that Z? is a diagonal matrix with non-negative diagonal coefficients and such that 
det(D) = 0. Let D' be a diagonal matrix with positive diagonal coefficients, and such that D ~ DD' and 
detD' = 2. We obtain 

Ki{tt oD) = Klin o {DD')) = 2'^Ki{tt o D), 
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which imphes that Ki{tt o D) = Q and concludes the proof. 



o 



The next proposition shows that the exponent p used for measuring the approximation error plays a 
rather minor role. By Kp we denote the error function associated with the exponent p. 

Proposition 2.2 There exists a constant c > such that for all 1 < pi < p2 < oo we have on 

cK < K < K < K 

Proof: For any function f E V = {l"^) and for any 1 < Pi < |52 < co by a standard convexity argument 
we obtain that 

||/|lLl(I<i) < ll/IUpjCI'') < ll/llip^pd) < ||/|U^(Id)- 

Using it follows that 

Ki < Kp^ < Kp^ < Xoo 
on Hm. Furthermore, the following semi norms on Hm 

kli Ik - I^IUiP'i) and |7r|oo := Ik - iTrH^^pd) 

vanish precisely on the same subspace of H^, namely Vi D Hm = {tt S : tt = Itt}. Since 
has finite dimension, it follows that they are equivalent. Hence, there exists a constant c > such that 
c| • loo < I • |i on Using it follows that cK^o < Ki^ which concludes the proof. o 



2.2 Desirable properties of the projection operator 



The examples of projection operators presented in §1.4| share some important properties which allow to 
obtain the explicit expression of the error function Kj. These properties are defined below and called 
Ha, or H^^,. They are satisfied when operator I is the interpolation at equispaced points (Definition 
1.12), at Tchebychev points (Definition |1.13[ ), and usually on the most interesting sets of ot her p oints. 



They are also satisfied when I is the L2(I^ orthogonal projection onto or (Definition 1.2). 

The first property reflects the fact that a coordinate Xi on I'^ can be changed to —Xi, independently 
of the projection process. 

Definition 2.3 {H± hypothesis) We say that the interpolation operator I satisfies the H± hypothesis 
if for any diagonal matrix D with entries in ±1 we have for all f £ V 

l{foD) = l{f)oD. 

The next property implies that the different coordinates xi, • • • ,Xd onY^ play symmetrical roles with 
respect to the projection operator. 

Definition 2.4 (iJ^- hypothesis) If is a permutation matrix, i.e. {Ma)ij := Sia{j) for some per- 
mutation a o/ {1, • • • ,d}, then for all f €z V 

I(/oM,) = I(/)oM,. 



According to (1.8), the projection operator I reproduces the space of polynomials P^,. However, in 
many situations the space Vi of functions reproduced by I is larger than P^.. In particular Vi = P" when 
I is the interpolation on equispaced or Tchebychev points, and Vi = Pfc (resp P^, Pj,*) when I is the 
^2(1'') orthogonal projection onto P/c (resp P^, P^^*). 

It is particularly useful to know whether the projection operator I reproduces the elements of P^, and 
we therefore give a name to this property. Note that it clearly does not hold for the ^2(1'') orthogonal 
projection onto P^. 

Definition 2.5 (iJ* hypothesis) The following inclusion holds : 

Pt c Vi. 
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On the contrary it is useful to know that some polynomials, and in particular pure powers s™, are 
not reproduced by I. 

Definition 2.6 (i?** iiypothesis) 

If J2 ^^^re'l^i then (Ai,--- ,Ad) = (O,--- ,0). 

l<j<d 

This condition obviously holds if !(/) G (polynomials of degree < A: in each variable) for all /. Hence, 
it holds for all the examples of projection operators given in the previous subsection §1.4| 

2.3 Explicit formulas 

In this section we provide the explicit expression for K when some of the hypotheses H±, H^, H^, or iJ** 
hold. Let TT € and let be the corresponding coefhcient of X™ in tt, for all 1 < z < d. We define 



i<i<fl 



and 



s(7r) := #{1 < i < d : A, > 0}. 



If ^^'^ is identified by (|1.16|) to an element of IEI,„, then one has 



ml 



m 



\ 



11 a^(-) 



Proposition 2.7 If m is odd and if H±, and iJ* hold, then 

Kp{TT) = C{p)K,{n), 

where 



C{p) := 



l<i<d 



, l<i<d 



> 0. 



Proposition 2.8 If m is even and if H^, iJ* and iJ** hold then 

Kp{n)=C{p,s{^))K,{n). 

Furthermore, 



C{p,0)^Cip,d) 



E ^r-i E 



, l<i<!i 



> 0. 



Other constants C{p, s) are positive and obey C{p, s) — C'{p, d — s). 



(2.20) 



(2.21) 



Next we turn to the proofs of Propositions |2.7| and |2.8| 
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Proof of Proposition 2.7 Let tt e E„i and let Xi be the coefficient of X™ in tt. Denote by 



l<i<d 



SO ttiat TT — TT* € and, more generally, tt o D — tt^, o D e for any diagonal matrix D. The hypothesis 
-ff* states that the projection operator I reproduces the elements of P^, and therefore 

TT o D - I(7r o Z?) = vr, o £) - I(7r* o D). 



Hence, Ki{'k) = Ki{n^,) according to ( |1.11[ ). If there exists ioi 1 < *o < rf, such that Xi^ — 0, then we 
denote by D the diagonal matrix of entries Dn — 1 if i iQ and if i — ig. Applying Proposition |2 . 1 1 we 
find 

Ki{tt) = Klin,) = Klin, o D) ^ [det D)"^ Kiin,) = 0. 

which concludes the proof. We now assume that all the co effic ients A^, I < i < d, are different from 0, 
and we denote by Ej be the sign of A^. Applying Proposition 2.1 to the diagonal matrix D of entries |Ai| 
we find that 



Kj{tt) = Ki{tt,) ^ (detL>)^if/(7r, oD-i) = K,{n)Ki ^ 



Using the H± hypothesis with the diagonal matrix D of entries Da = Si, and recalling that m is odd, we 
find that 



We now define the functions 
It follows from ( |1.11| that 



\l<i<d J \l<i<d 

- I(Xr) for 1 < i < d. 



Kr I E 

, l<i<d 



inf 



l<i<d 



where the infimum is taken over all c?- vectors of positive reals of product 1. Let us consider such a d- vector 
(ai, • • • , ad), and a permutation a of the set {1, ■ ■ • ,d}. The i/o- hypothesis implies that the quantity 



l<i<d 



is independent of a. Hence, summing over all permutations, we obtain 



-y 



1 

> - 
- d 



E 3^ 

l<i<d 



E 



(2.22) 



The right-hand side is minimal when ai = • • • = = 1, which shows that 



E ^i9i 


> 


E 9^ 


= C{p) 


l<i<d 




l<i<d 





with equality when = 1 for all i. Note as a corollary that 

Klin,) = IItt^ - I(7re)||i^(id) = C{p) where tt^ = E ^^^^ 



(2.23) 



l<j<d 



11 



It remains to prove that C{p) > 0. Using the hypothesis H±, we find that for all <E {±1} we have 



l<i<d 



dp). 



In particular, for any \ <iQ < d one has 



2|l5»olUp(i'') < 




+ 




< 2C{p) 




l<i<d 




l<i<d 





If C{p) = 0, it follows that g,„ = and therefore that X:;^ = HXl'^), for any I < io < d. Using the 
assumption H^,, we find that the projection operator I reproduces all the polynomials of degree m = k + 1, 



which contradicts the definition (1.8) of the integer k 



Proof of proposition 2.8 We define A^, tt* and e.i G {±1} as before and we find, using similar 
reasoning, that 



, i<j<£; 



For 1 < s < d we define 



s+l<i<d 



From the hypothesis it follows that Ki{tt) — K^,{Tr)C{p, s{tt)). 

Using again H^^ and the fact that Kj{tt) = Kj{—Tr) for all tt e H™, we find that 



K, 



J2 xr- E 



C{p,d-s). 



, l<i<s 



s+l<i<a 



, l<i<d-s 



-s + l<i<d 



We define gi := X™ — I(X™), as in the proof of Proposition 2.7 We obtain the expression for C(p, 0) 



by summing over all permutations as in (2.22) 



C(p,0) 



E^ 

l<i<d 



This concludes the proof of the first part of Proposition 2.8 We now prove that C{p,s) > for all 



1 < p < oo and all s G {0, • • • , d}. To this end we define the following quantity on 



Hk 



E 



l<i<d 



Note that \\a\\K — if and only if 



E «.xr= E 



l<i<d 



l<i<d 



and the hypothesis precisely states that this equality occurs if and only if — 0, for all 1 < i < d. 
Hence, || • \\k is a norm on M''. Furthermore, let 



n 



l<i<d 
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Then 

C(p, s) = inf \\a\\K- 

a^Es 

Since Eg is a closed subset of M'', which does not contain the origin, this infimum is attained. It fohows 
that C(p, s) > 0, and that there exists a rectangle of unit volume such that 

Kii^e) - he ~ lT^e\\L,(R,) = C[p,s{n,)) where TT, = ^ £,Xr. (2.24) 

l<i<d 

O 



3 Proof of the approximation results 

In this section, let the block i?o, the integer m, the function / G C""(i?o) and the exponent p be fixed. 
We conduct our proofs for 1 < p < oo and provide comments on how to adjust our arguments for the 
case p = oo. 

For each a; G Mq by /i^ we denote the m-th degree Taylor polynomial of / at the point x 

= M.(x) E a^(-) ^^r^' (3-25) 

|a| <m 



and we define tTx to be the homogeneous component of degree m in fi^ 

ax"^""^ a! 



-^ = -^w- E 7d^-)^- (3-26) 

|a|— m 

Since tt^: and fx^ are polynomials of degree m, their m-th derivative is constant, and clearly d"^T:.j. — 
= (i™/(a;). In particular, for any x G Rq the polynomial /x^. — tt^; belongs to (recall that 
k = m — I) and is therefore reproduced by the projection operator I. It follows that for any x d Rq and 
any block R 

T^x —^rt^x = ^J-x —^RfJ-x- (3.27) 
In addition, we introduce a measure p of the degeneracy of a block R 

, , diam(i?)'' 



Given any function g S C™'{R) and any x € R we can define, similarly to (3.26), a polynomial Tr^ G 
associated to g at x. We then define 

Wd^'ghuR) sup ( sup |7r,(u)| ) . (3.28) 

x€R y|M| = l / 

Proposition 3.1 There exists a constant C ~ C{m, d) > such that for any block R and any function 
g G C""(i?) 

\\9 ^lR9\\L,iR) < C\R\^ piR)'^\\d"'9\\L^iR)- (3.29) 

Proof: Let xq ^ R and let go be the Taylor polynomial for g of degree m — 1 at point xq which is defined 
as follows 



9oiX):= E 



d^fixo) (X - xoT 



dx°' a\ 

\a\'Cm—l 



Let x £ R and let x{t) = xo + t{x — xq). We have 

nl jm. 



I \ ( \ ^ [ dJ^9x{t){x-xa) 
g{x)^go{x)+ — (1-t) 



„, dt. 

t=o 
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Hence, 



\g{x)-go{x)\ < I ||d™5llL^(fl)k-xor"(l-t)™dt< ^^||d"g||L^(fl)diam(i?r 



(3.30) 



Since go is a polynomial of degree at most m — 1, we have g^ — IgQ. Hence, 

ll3-ii?.g||Lp(fl) < \R\h\g~'iR9\\L^{R^) 

= l^l'^ll(5"5o) -ifl(.g-.go)llL^(i?.) 

< (l + Ci)|i?|?||.g-.go||L^(fl), 



where Ci is the operator norm of I : F — > Combining this estimate with (3.30), we obtain (3.29). 



3.1 Proof of Theorem 1.4 (Lower bound) 



The following lemma allows us to bound the interpolation error of / on the block R from below. 
Lemma 3.2 For any block R d Rq and x £ R we have 

\\f-Uf\\L,(R) > \R\^ (i^/(^.)-c^(diami?)p(i?)^), 
where the function lo is positive, depends only on f and m, and satisfies a;((5) — > as i5 — > 0. 



Proof: Let h :— f — ii^, where fi^ is defined in (3.25) Using (3.27), we obtain 



/- Ifl. /II Lp(fl) > Ikx - Ifl. 7r2:||Lp(_R) - 11^ - li?, ft-||Lp(_R) 

> \R\^Kii7r,)^\\h-lRh\\L^(^n^, 



and according to (3.29) we have 

\\h~lnh\\L^^i,)<Co\R\hiR)'^\\d"'hU^^j,y 

Observe that 

IM'^/^IIl^W = ||d™/-d'"^.||L^(i?) = |M'"/-d™/(a:)|U^(«). 
We introduce the modulus of continuity UJ.^, of the m-th derivatives of /. 



w*(r):= sup ||d'"/(a;i)-d"7(2:2)|| - sup sup |7r,, (u) - (u)| (3.31) 

xi,X2&Ro- xi,X2£Ro- y|M|<l J 

\xi—X2\<r \xi—X2\<r 

By setting w = Co we conclude the proof of this lemma. o 

We now consider an admissible sequence of block partitions {TIn)n>o- For all > 0, i? € TZn and 
X G R, we define 

(I)n{x) := |i?| and i/^Nix) := {Ki{Tr^) - a;(diam(i?))p(i?)T)^ , 

where A+ := max{A,0}. We now apply Holder's inequality /^^ /1/2 < ||/i||Lpjj?o)ll/2|lLp,(i?o) with the 
functions 



fi^(l>N and /; 



and the exponents pi = ^ and p2 = ^ . Note that ^ 



2 = (Pn 



Pl P2 \p d 



1. Hence, 



Rq 



Ro 



(3.32) 
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Note that fj^^ (jy^^ — #(7?. at) < N. Furthermore, if R £TZn and x £ R then accordmg to Lemma 

(t>N{x)'^^N{x) = \R\^-hNix) < \Rr--\\f'lRf\\L,(R)- 

Hence, 



3.2 



< 



- 11/ - I-R /IILp(_Ro)- 



Inequahty (3.32) therefore leads to 



p(flo) 



(3.33) 



(3.34) 



Since the sequence {TZn)n>o is admissible, there exists a constant Ca > such that for all N and all 
i? e TZn we have diam(i?) < CaN~^. We introduce a subset of 72.^ C TZn which collects the most 
degenerate blocks 

n'j, = {RenN: piR) > uj{CaN--^)-^^}, 



where ui is the function defined in Lemma 3.2 By i?'^ we denote the portion of Rq covered by TZ'^^. For 
all X G Ro\ R'j^ we obtain 

i^N{x) > Ki{tt,)-lo{CaN-"^Y-"^. 
We define Sn '■= uj[CaN^^Y^^ and we notice that eat — > as iV — > oo. Hence, 



||V'Jv|II^(B„) > ||(i^/(vr:,) -eAr)+||^ 



> (if/(7r,)-£Ar), li: .~G-\R' 



LARo\R'^) 

T 



where C :— maxx^Rg Kj{7Tx). Next we observe that \R']^\ — > as iV — +oo: indeed for all R E TZ'p^ we 
have 

\R\ = diam{Rfp{R)~^ < CiN~^Lo{CAN-i)^^ . 

Since #(7?.'^) < A'', we obtain \R'j^\ < Cj^uj{CAN^i)'^ , and the right-hand side tends to as A'' — oo. 
We thus obtain 

ljm|nf ||V'w||L,(flo) >^lini^||(^i'/(7ra;)-£Ar)+||^^(^^) = \\Ki{tTx)\\lar„). 



Combining this result with (3.34 1, we conclude the proof of the announced estimate. 
Note that this proof also works with the exponent p = cxd by changing 

mp \ P m 

<t^N into UA^n\\1^(b,„) 



in (3.321 and performing the standard modification in (3.331. 



Lp(Ro,n) 



Remark 3.3 As announced in Remark 1.6 this proof can he adapted to the weighted norm 
associated to a positive weight function \ l G C'^{Rq) and defined in ( 1.1!^ . For that purpose let rjv ■— 
sup{diam(i?) : R e TZn} and let 

Qn{x) := inf D,{x'). 

x'eRa 

The sequence of functions fl n increases with N and tends uniformly to as N ^ oo. If R E TZn o,nd 
X G R, then 

11/ - Ii? /llLp(fl,n) > ^N{x)\\f - iRfh^R)- 

The main change in the proof is that the function ippf should be replaced with -tj/j^ := fijvV'iV- Other details 
are left to the reader. 
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3.2 Proof of the upper estimates 



The proof of Theorems 1.5 (and 1.7 1 is based on the actual construction of an asymptotically optimal 



sequence of block partitions. To that end we introduce the notion of a local block specification. 

Definition 3.4 (local block specification) A local block specification on a block Rq is a (possibly 
discontinuous) map x i-> R{x) which associates to each point x € Rq a block R{x), and such that 

• The volume \R{x)\ is a positive continuous function of the variable x ^ Rq. 

• The diameter is bounded : sup{diam(i?(a;)) : x G Rq} < oo. 

The following lemma shows that it is possible to build sequences of block partitions of Rq adapted in 
a certain sense to a local block specification. 

Lemma 3.5 Let Rq be a block in J?*^ and let x i— > R(x) be a local block specification on Rq. Then there 
exists a sequence {Vn)n>i of block partitions of Rq, Vn = U V^, satisfying the following properties. 

• (The number of blocks in Vn is asymptotically controlled) 

lim = / \R{x)\-^dx. (3.35) 

• (The elements ofV^ follow the block specifications) For each R G there exists y £ Rq such that 

R is a translate ofn^^R{y), and \x — y\ < — ^ — — for all x (£ R. (3.36) 



(The elements ofV^ have a small diameter) 



lim sup diam(i?) = 0. (3.37) 

n-i-oo \^ ^gp2 j 

Proof: See Appendix. o 
We recall that the block i?o, the exponent p and the function / € C™{Rq) are fixed, and that at each 



point X £ Rq the polynomial G H,„ is defined by (3.26). The sequence of block partitions described 



in the previous lemma is now used to obtain an asymptotical error estimate. 
Lemma 3.6 Let x i~> R{x) be a local block specification such that for all x £ Rq 

Wt^x - ^r{x){'^x)\\l^(r(x)) < 1- (3.38) 
Let {Vn)n>i be a sequence of block partitions satisfying the properties of Lemma 3.5, and let for all N > 

n{N) := max{n > 1 : #('P„) < N}. 
Then TZn :— 'Pn(N) *s an admissible sequence of block partitions and 

limsupiV?||/-l7j„/|U^(fl„) < f / R{x)-Hx\ . (3.39) 

Af->-oo \J Ra / 



Proof: Let n > and let R^Vn- If i? e V\ then let y e i?o be as in ( |3.36[ ). Using ( |3.29[ ) we find 

ll/-Ii?/IUp(ii) < Iky - + IK/ - TTy) - Ifl(/ - 7ry)||i^(fl,) 

< Ik, - lR(y) ^y\\L,(Riy)) + C\R\^ diam(i?)™||d™/ - d^iTyU^^^^ 

< n-^ +Cn-^|i?(y)|^diam(i?(y))"||d"7-d"7(y)l|L^(fl.) 

< n-^(l + C"a;*(n"Miani(i?o))), 
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where we defined C := Csupygg^ \R{y)\p diam(i?(y))™, which is finite by Definition 3.4 We denoted by 
oj* the modulus of continuity of the m-th derivatives of / which is defined at (3.31 ). We now define for 
all n > 1, 

6n '■= sup diam(_R). 



According to (3.371 one has (5„ ^ a s n oo. If i? G T'^, then diam(i?) < n and therefore 
\R\ < diam(i?)'' < n^^'^S^. Using again (|3.29|), and recalling that 7 = + p ^^^^ 



ll/-Ii^/llL,(i^) < diam(i?)"||d"/|U^,,^, < C"n-^<5„^ 

where C" = C||d™/||L^(^g) ■ From the previous observations it follows that 

ll/-Ip„ /IU,(i^o) < max Wf-lRfhAR) < #{Vn)^n-'^ max{l + C'c^, (n"! diam(i?o)), C's!}. 

Hence, 

limsup#(7'„)-^n^||/-Ip„/|U^(fl„) < 1. 
Combining the last equation with ( 3.35[ ), we obtain 



limsup#(7'„)^||/-Ip„/|U^(Ko) < / R{x)-^dx 



Ro 



The sequence of block partitions TZn :— 'Pn{N) clearly satisfies ^{TZn)/N — 1 as — cx) and therefore 
leads to the announced equation (3.391. Furthermore, it follows from the boundedness of diam(i?(a;)) on 



i?o and the properties of Vn described in Lemma 3.5 that 



sup ( #{Vn)^ sup diam(i?) I < 00 
n>i \ Rev„ J 

which implies that 7?. at is an admissible sequence of partitions. o 
We no w ch oos e ad equate local block specifications in order to obtain the estimates announced in 



Theorems 



1.5 



and 



1.7 



For any M > diam(I'') ^ y/d wc define the modified error function 
Km{tt):^ inf ||7r-Ifl7r||i^(fl), 

|K| — 1, 
diam(_R)<M 



(3.40) 



where the infimum is taken on blocks of unit volume and diameter smaller that M. It follows from 
a compactness argument that this infimum is attained and that Km is a continuous function on Hm. 
Furthermore, for all tt G Hm, M Km{t^) is a decreasing function of M which tends to Kj{tt) as 
M 00. 

For all X Cz Rq we denote by R*j^j{x) a block which realises the infimum in i^jvf (tTj.). Hence, 
\R*Mix)\ = 1, dia.m{Rlj{x)) < M, and Km{iTx) = hx ~ ^r.Ij(x) t^x\\lp(ri,(x)) 
We define a local block specification on i?o as follows 

Rm{x) {Km{^x) + M-^^Rljix). (3.41) 

We now observe that 

\\t^x - ^Rm{x) T^x\\LpiRM{x)) = Km{t^x){Km{t^x) + Af"^)"^ < 1. 



Hence, according to Lemma 3.6 there exists a sequence {TI^)n>i of block partitions of Rq such that 



limsup7V7||/-l7,M/|U^(;j„) < \\KMi7rx)+M-'\\LARo)- 

N—>-oo 
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Using our previous observations on the function Km, we see that 



lim \\Km{^.) + M-^\\L^(n.o) - \\Ki{tiMlaRo)- 

M —¥oo 

Hence, given £ > we can choose M {e) large enough in such a way that 

l|i^Af(e)(7r.) + M(£)-'||L.(%) < ||i^/(^.)|U,(fl„)+£, 



which concludes the proof of the estimate (1.18) of Theorem 1.5 



For each N let M = M{N) be such that 

^^ll/-lK-/IU,(flo) < II^m(^x) + M-1|U^(j,„)+M-i 
and M{N) oo as N — ^ oo. Then the (perhaps non admissible) sequence of block partitions TZn :- 



n 



M(N) 
N 



satisfies (1.171 which concludes the proof of Theorem 



1.5 



We now turn to the proof of Theorem |1.7[ which follows the same scheme for the most. There exists 
d functions Xi{x), • • • , Arf(x) G C'^{Ro), and a function x i— >■ tt^{x) S such that for all x G Rq we have 

l<i<d 

The hypotheses of Theorem |l.7| state that Ki (^^;if ) = Ki{'Kx) does not vanish on Rq. It follows from 
Propositions 2.7 and |2.8| that the product \i{x) ■ ■ ■ Xd{x) is nonzero for all x G Rq. We denote by 



Si e {±1} the sign of A^, which is therefore constant over the block Rq, and we define 



l<i<d 



The proofs of Propositions 2.8 and 2.7 show that there exists a block i?^, satisfying = 1, and such 
that Kj^TT^) — Wtt — Iji^ 7r||i (^^j. By D{x) we denote the diagonal matrix of entries |Ai(a;)|, • • • , |A(i(a;)|, 
and we define 

R*{x) := {detDix))^D{x)-^R^. 



Clearly, |i?*(x)| — 1. Using (1.5) and the homogeneity of tt^ G Ttlm, we find that 

hx-'^R'ix) T^x\\LpiR-{x)) = idetD{x))iKj{T:i;) ^ Ki{tt^). 
We then define the local block specification 

R{x) ■.= Ki{tt,)-^^R*{x). 



(3.42) 



The admissible sequence {TZn)n>i of block partitions constructed in Lemma 3.6 then satisfies the optimal 
upper estimate (1.17), which concludes the proof of Theorem 1.7 o 



Remark 3.7 (Adaptation to v^reighted norms) Lemma 3.6 also holds if (3.38) is replaced with 

Vl{x)\\tt^ - Ir{x){t^x)\\Lp{R{x)) < 1 

and if the Lp(i?o) norm is replaced with the weighted Lp{Ro,fl) norm in ( 3.3£^ . Replacing the block 
Rm{x) defined in (3.4-1) with 

r'm{x) -.^nixy^RMix), 



on e can easily obtain the extension of Theorem 1.5 to weighted norms. Sim ilarly , replacing R{x) defined 
in (3.4-2) with R'ix) :— VL{x)~ R{x) , one obtains the extension of Theorem 1.1 to weighted norms. 

APPENDIX 
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Figure 1: (Left) the initial uniform (coarse) tiling Q3 of i?o- (Right) the set of blocks in green and 
the set of blocks V^* in red. 



A Proof of Lemma 13.5 



By Qn we denote the standard partition of Rq e IR*^ in n'^ identical blocks of diameter n~^diam(i?o) 
illustrated on the left in Figure 1. For each Q e Q„ by xq we denote the barycenter of Q and we consider 
the tiling Tq of M'^ formed with the block n~^R{xQ) and its translates. We define V}^{Q) and as 
follows 

VlM ■.= {R^Tq: RdQ} and Vl := |J Vl{Q)- 

QeQ„ 

Comparing the areas, we obtain 



\Q\ 



2d 



QeQ„ 



QeQ, 



\n-^R{xQ)\ 



E \Q\\Ri^Q)\-'- 

QeQ„ 



From this point, using the continuity of a; h- !■ |i?(a;)|, one can easily show that 



J^jR{x)\-'dx 

as n — > cxD. Furthermore, the property (3.36) clearly holds. In order to construct V^, we first define two 



U ^niQ)- 



sets of blocks V^* (Q) and V^* as follows 

7?2*(Q) {i?n Q : i? e Tq and i? n 9Q 7^ 0} and V^* 

Comparing the surface of dQ with the dimensions of R{xq), we find that 

#{V?l*{Q))<Cn''-' 

where C is independent of n and of Q e Q„. Therefore, 4h{Vf*) < Cri^'^^^. The set of blocks V^^ is then 
obtained by subdividing each block of Vf* into o{n) (f or ins tance, [ln(n)J identical sub-blocks, in such 
a way that #(7^^) is o{n^'^) and that the requirement (3.37) is met. 
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